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A  theoretical  analysis  is  presented  of  the  heterogeneous  opposed 
flow  diffusion  flame  (HOFD) .  In  such  a  system  a  gaseous  reactant  (G) 
impinges  on  a  solid  reactant  (S)  undergoing  sublimation  (or  vaporiza¬ 
tion)  and  a  diffusion  flame  is  established  at  the  stagnation  point. 

With  increasing  mass  flux  of  the  gaseous  reactant  (m  )  the  mass  flux  of 

G 

solid  reactant  (m  )  grows  until  the  total  consumption  rate  of  reactants 
b 

reaches  a  limiting  value  due  to  limitations  imposed  by  the  kinetics  of 
reaction.  The  analysis  of  the  axi-symmetrical  stagnation  point  flow  with 
chemical  reaction  allows  evaluation  of  the  kinetics  of  chemical  reaction 
at  temperatures  of  interest  to  combustion  processes.  The  theoretical 
development  is  applied  to  the  solid-propellant  deflagration  of  an  ammon¬ 
ium  perchlorate-catalyst-fuel  system  and  the  pyrolysis  of  solid  polymers 
in  an  oxidizing  atmosphere. 
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Introduction 


For  the  study  of  reaction  kinetics  of  gas  mixtures  at  tempera¬ 
tures  of  interest  to  combustion,  the  opposed  flow  diffusion  flame  (OFD) 
has  been  particularly  valuable.  In  experiments  employing  OFD,  one  mea¬ 
sures  the  highest  mass  flux  rates  of  opposing  gas  jets  of  fuel  and  oxidizer 
that  can  be  consumed  in  a  diffusion  flame  located  in  the  stagnation  re¬ 
gion.1  A  series  of  analytical  studies2-6  examined  the  relationship  be¬ 
tween  the  condition  of  diffusion  flame  extinction,  the  maximum  kinetics 
of  gas-phase  reaction,  and  the  flame  propagation  velocity  of  premixed 
gases.  These  theoretical  developments  were  concerned  with  homogeneous 
OFD  in  which  the  two  reactants  were  introduced  at  the  same  mass  flow 
rate  per  unit  area  into  the  reaction  region.  To  apply  this  technique 
to  solid-propellant  combustion  kinetics  we  have  treated  a  theoretical 
model,  the  heterogeneous  opposed  flow  diffusion  flame  (HOFD) ,  in  which 
a  gaseous  fuel  jet  impinges  on  a  solid  oxidizer  surface,  such  as  ammon¬ 
ium  perchlorate  (AP) ,  which  undergoes  sublimation.  This  model  considers 
exothermic  chemical  reaction  in  the  diffusion  flame  and  heat  losses  to 
the  environment . 

The  fundamental  difference  from  the  homogeneous  OFD  is  the  con¬ 
dition  at  the  solid/gas  interface,  the  source  of  the  oxidizer  species 
in  the  case  of  a  solid  propellant  based  on  AP.  Its  steady-state  surface 
properties  are  no  longer  an  independent  variable;  rather  the  mass  flux 
of  reactant  leaving  the  solid  surface  is  a  function  of  the  heat  and  mass 
transport,  and  the  thermodynamic  properties  of  the  solid.  Altnough 
similar  in  some  respects  to  the  opposed-flow  model  chosen  by  Spalding2, 
our  analysis  of  the  HOFD  postulates  a  much  more  realistic  velocity  field 
which  satisfies  the  conditions  at  upstream  infinity. 
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Mathematical  Formulation 


As  shown  schematically  in  Figure  1,  the  analysis  concerns  the 
opposed  jet  system  and  involves  a  stream  of  fluid  (the  gaseous  fuel, 
subscript  G)  moving  in  the  negative  z  direction  toward  a  perpendicular 
flat  surface  comprising  the  solid  reactive  oxidizer  (subscript  S)  under¬ 
going  sublimation.  The  axi-symmetrical  stagnation-point  flow  leads  to 
the  following  equation  for  conservation  of  energy  (with  the  velocity 
components  defined  in  Figure  l) : 

C  [(puT)  +  puT/r  +  (pvT)  ]  =  X  (T  +  T  /r  +  T  )  +  QR  (l) 

G  r  z  rr  r  zz  P 

where  C  is  the  specific  heat  of  the  gas  at  constant  pressure,  p  its 
G 

density,  X  its  thermal  conductivity,  Q  the  exothermic  heat  of  reaction, 

and  Rp  the  chemical  rate  of  product  generation.  Throughout  the  analysis, 

the  Lewis  number  is  set  equal  to  unity.  ( X/c  pD  =  l) ,  and  the  variable 

G 

subscripts  r  and  z  indicate  partial  differentiation. 

The  equation  for  conservation  of  total  mass  is 

(pu)r  +  pu/r  +  (pv)z  =  0  ,  (2) 

while  that  for  the  mass  of  species  K  reads 


(puVK  '  pDVr  *  <puVPDYKr)/r 


*  (PVYK  -  PDVz  *  RK 


(3) 


where  Y  is  the  mass  fraction  of  species  K  and  R  is  its  rate  of  forma- 
K  K 

tion  in  the  reaction.  The  species  are  taken  to  be  oxidizer  X,  fuel  F, 

product  P,  and  inert  gas  I.  By  introducing  the  constants  w  =  R  /QR 

K  K  P 

we  can  simplify  the  equation  to 


p“  <YK-»KT)r  +  pv  <YK-"KT)z 


=  PD  [(YK  -  *K  T)rr  +  (YK  '  »K  TVr  *  (YK  ‘  VW  '  (4) 
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reads 


(9) 


\  T  (y')2  +  Q  R  -  0 

yy  p 

Solution  of  the  Equations 

It  is  now  assumed  that  (y')2  has  a  reasonable  form  as  a  function 
of  z  so  that  the  distribution  of  reaction  rate  can  take  that  form: 

Rp  =  (*j/q)  (y')2  (10) 

where  J  is  a  constant.  Then 

T  =  A  +  By  -  Jy2/2  (11) 

where  A  and  B  are  constants.  Thus  all  dependent  variables  have  explicit 
forms  as  functions  of  y  on  r  =  0. 

A  reasonable  form  for  y  is 

y  =  tanh  az  (12) 

with  some  constant  a.  One  finds  on  r  =  0  by  Eq .  (7) 

pv  =  -2apD  tanh  az  (13) 

so  lhat  to  fit  the  condition  at  G 

a  =  m  /2pD  (14) 

u 

An  extension  of  pv  to  points  off  the  axis  of  symmetry  which  gives 
a  jet-like  flow  is 

pv  =  -  m^  tanh  az  sech2(br2/2)  (15) 

where  constant  b  can  be  related  to  the  jet  width. 

The  equation  for  conservation  of  total  mass  can  vhen  be  solved 

to  give 

pu  =  (am  /br)  sech2  az  tanh  (br2/2)  (16) 

G 
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This  solution  satisfies  the  boundary  condition  at  G,  pu  =  0  at  z  =  °=, 

but  at  the  solid  surface  (z  =  z  ) ,  pu  =  0  only  on  r  =  0.  Since  the 

s 

solution  is  needed  on  r  =  0,  this  error  is  not  likely  to  affect  the  re¬ 
sults  seriously. 

The  solution  for  T  on  r  =  0  satisfying  the  condition  at  G  is 


T  =  Tg  +  B*(l-y)  -  J(l-y)2/2  (17) 

where  B'”'  is  a  convenient  form  for  the  remaining  constant  of  integration. 
The  conditions  at  the  solid  surface  determine  y  =  -m  and 

{ ( 1+m)  2  J  [m  Cg+(l-m)CG]  -  2m  [L  +  C^.-T^]} 

B  (l+m)  [2m  C  +  (l-m)  C  ] 

S  G 

where  m  =  m  /m .  This  ratio  is  restricted  to  0  <  m  <  1  by  Eq .  (15). 

b  VJ 

Boundary  conditions  at  S  and  G  similarly  determine  the  constants 

of  integration  A  and  B  in  the  solution  forms  for  the  Y  . 

K  K  K 

The  maximum  temperature  is 

T  =  T  +  B-"-2/2J  (19) 

MG7  ' 

and  it  occurs  at  y  =  y  .  where  y  =  1  -  8"”'/j.  It  can  be  readily  deter- 
mined  from  the  reaction  rate  R  at  point  M: 


J  =  4\Q  R  /m2  (1 
PM7  G 


-  y2)2 
■V 


ilication  to  Heterogeneous  Combustion 


For  the  special  case  for  which  (l)  the  reaction  rate  is  of  first 
order  in  the  fuel  and  oxidizer  concentrations,  (2)  the  conductive  heat 
loss  into  the  solid  phase  [c  (T  -T  )]  is  negligible  relative  to  the 

o  o  A 

heat  of  sublimation  L,  and  (3)  the  temperature  maximum  is  located  at 
the  stagnation  point  (y  =  0) ,  the  equations  just  derived  become  greatly 
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simplified.  For  such  a  case  the  rate  of  product  generation  is  given 

by 


Rp  =  P2  (Yp)  (Yx)  Z  exp  (-E/RT) 


(21) 


where  Y  and  Y  are  the  weight  fractions  of  fuel  and  oxidizer,  E  the 
F  X 

activation  energy ,  R  the  gas  constant ,  and  Z  the  preexponential  term 

in  the  Arrhenius  expression.  This  equation  applies  equally  at  the 

stagnation  point  where  R  -*  R  and  T  -*  T  .  To  simplify  the  calcula- 

tion,  we  make  the  additional  approximation  that  at  this  point  Y  and 

F 

Y  are  unity.  Under  these  conditions  Eq.  (20)  reduces  to 
X 


=  4p  2AQZ  exp  (-E/RTM)/m2G 


(22) 


Finally,  in  combination  with  Eqs.  (18)  and  (19),  after  neglecting  mCc 

with  respect  to  (l-m)C  ,  one  obtains 

G 


m-  =  202XQ(Cr_/L)  (1-m2)  Z  exp  ("E/rT^ 


'g 


(23) 


where  T  =  T  +  [l/c  (l-m2)] 

MG  G 

The  function  given  by  Eq.  (23)  has  the  properties  depicted  in 

Figure  2.  With  increasing  mass  flux  of  the  fuel  component  (m  ) ,  the 

G 

steady-state  mass  flux  of  solid  oxidizer  (m  )  grows  until  it  attains  a 

b 

maximum  value  at  the  limiting  ratio  m*  =  m*/m*.  Further  increase  in 

S  G 

m^  causes  the  kinetics  of  reaction  to  limit  further  consumption  of  re¬ 
actants  and  leads  to  extinction.  At  the  maximum  point  d  m  /d  m  =0 

S  G 

or  d  m  /  d  m  =  -m  /m  so  that 
G  G 


E7  = 


2m’:' : 


(1+I/-I#2) : 


m 


f;:' 2  ( 1  -i#  2 ) 


(24) 


where  E7  =  E/RT  and  L7  =  l/cT  . 

G  G  G 


In  Table  1  we  have  computed  E7  for  different  values  of  m  and  L7. 


Also  for  values  of  m*  2  near  unity  Eq.  (24)  reduces  to 

E '  Sr  l'/ ( 1-m'*2)  (25) 

It  is  apparent  that  the  activation  energy  for  the  fuel-oxidant 

gas-phase  reaction  may  be  evaluated  from  the  critical  mass  flux  ratio 

m*.  An  interesting  application  of  this  theoretical  model  is  to  be 

found  in  the  pyrolysis  measurements  of  polymethyl  methacrylate  (PMM)  by 

McAlevy  and  coworkers.7  In  an  experimental  study  involving  a  stream  of 

oxygen  impinging  on  the  solid  polymer,  the  linear  regression  rate  of 

the  polymer  attained  a  maximum  corresponding  to  a  critical  mass  flux 

ratio  of  m'"'  =  m'::'/m'"  =  1.5  x  10  2.  Based  on  an  effective  endothermic 
S'  G 

heat  of  sublimation7  of  about  1000  cal/g  one  calculates  by  substitution 
into  Eq.  (25)  an  activation  energy  of  8  kcal/mole  for  the  gas-phase 
combustion  process  of  PMM-oxygen.  As  shown  in  a  subsequent  publication8 
we  have  used  the  HOFD  theory  in  an  evaluation  of  the  reaction  kinetics 
of  the  system  ammonium  perchlorate-catalyst-hydrocarbon.  Further 
theoretical  work  is  in  progress  to  remove  the  restriction  of  i#  <  1 
from  the  analysis,  so  that  the  model  may  have  applicability  to  a  wider 
range  of  combustion  systems. 
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Table  1 


COMPUTATION  OF  E7  (Eq.  24) 


y/  L 

e' 

0.1 

0.3 

1 

3 

10 

0.6 

2.08 

1.35 

1.63 

3.23 

9.0 

0.8 

3.38 

3.11 

5.40 

12.8 

39.0 

0.9 

3.55 

4.71 

10.7 

29.3 

90.5 

0.95 

4.26 

7.70 

20.8 

58.8 

193.0 
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